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Abstract 



We investigate the evolution of different measures of "Gravitational En- 
tropy" in Bianchi type I and Lemaitre-Tolman universe models. 

A new quantity behaving in accordance with the second law of thermody- 
namics is introduced. We then go on and investigate whether a quantum cal- 
culation of initial conditions for the universe based upon the Wheeler-DeWitt 
equation supports Penrose's Weyl Curvature Conjecture, according to which 
the Ricci part of the curvature dominates over the Weyl part at the initial sin- 
gularity of the universe. The theory is applied to the Bianchi type I universe 
models with dust and a cosmological constant and to the Lemaitre-Tolman 
universe models. We investigate two different versions of the conjecture. First 
we investigate a local version which fails to support the conjecture. There- 
after we construct a non-local entity which shows more promising behaviour 
concerning the conjecture. 



I. INTRODUCTION 

The physics of the arrow of time in our universe seems somehow to involve gravitation 
|T|]. The origin of this arrow may be searched for in the initial state and early evolution of 
the universe 0]. 

The isotropy of the cosmic background radiation shows that the universe was to a very 
high degree in a state of thermodynamic equilibrium already 300.000 years after the Big 
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Bang. At the present time there are great thermal differences between matter in the stars 
and that in the interstellar clouds. Matter has developed away from a state of thermal 
equilibrium instead of towards it, in contradiction to what one could expect from the second 
law of thermodynamics. This is clearly due to gravity, which increases density fluctuations 
and hence increases the cosmic temperature differences. 

One would like to incorporate this tendency of gravity to produce inhomogeneities into 
a generalized second law of thermodynamics. Then one needs a quantity representing the 
entropy of a gravitational field. This quantity should vanish in the case of a homogeneous 
field and obtain a maximal value given by the Bekenstein-Hawking entropy of a black 
hole, for a gravitational equilibrium configuration in the form of the field of a black hole. 

One suggestion in this connection was Penrose's formulation of what is called the Weyl 
curvature conjecture{W CC) Wainwright and Anderson express this conjecture in 
terms of the ratio of the Weyl and the Ricci curvature invariants. 



(1) 



According to the conjecture vanishes at the initial singularity of the universe. The 
physical content of the conjecture is that the initial state of the universe is homogeneous 
and isotropic. As pointed out by Rothman and Anninos the entities and Cap-ysC"^"'^ 
are local entities in contrast to what we usually think of entropy. In this paper we will 
investigate various entities in the Bianchi type I and the Lemaitre-Tolman models. We will 
give a careful investigation of both the entities Cap-^sC"'^^^ and P^, and show that both 
of these fail to describe a behaviour according to the WCC We will therefore introduce a 
non-local entity which shows a more promising behaviour concerning the WCC. This entity 
is also constructed in terms of the Weyl tensor, and it has therefore a direct connection 
with the Weyl curvature tensor but in a non-local form. The local version of the WCC is a 
much stronger restriction on the cosmological models and it is not strange that it fails for a 
general model. 

The vanishing of the Weyl curvature tensor at the initial singularity is a very special initial 
condition. Due to the evolution during the inflationary era, however, this condition may 
be relaxed. A much larger variety of initial conditions are consistent with the homogeneity 
of the universe, as observed in the cosmic background radiation, in inflationary universe 
models than in models without inflation. As we shall see from the quantum calculations, 
an inhomogeneous universe is more likely to be spontaneously created than a homogeneous 
one. A universe with a large cosmological constant is also more likely to be created than a 
universe with a small cosmological constant. 

In sections II and III, respectively, we present a classical cosmological investigation of 
Ca/j-y^C"^^^ and P^ for the Bianchi type I and Lemaitre-Tolman universe models. Both of 
these entities fail to have the right entropic behavior. We therefore propose in section IV an 
entity which is non-local and shows a much better behaviour concerning the WCC. It should 
be noted however that this entity is by no means the true entropy for the gravitational field 
but it does have the right entropic behaviour. In section V we give a quantum cosmological 
treatment of the WCC in the stage set by the classical investigations. 
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II. THE WEYL CURVATURE CONJECTURE FOR A BIANCHI TYPE I 

MODEL: LOCAL VERSION 



Based on the results obtained in an earlier paper where the Bianchi type I minisuperspace 
model was investigated 0, we will consider the Weyl curvature conjecture (WCC) for a 
Bianchi type I model. Writing 

ds^ = -dt" + e^" [e^f] dx'dx^ (2) 

where (5 = diag(/3+ + v/3/5-,/3+ — -\/3/5_, —2/5+) the Weyl curvature invariant and the Ricci 
square turn out to be: 

Ca/sysC''^^' = ^1201 + + 3{(3l + i3l) + 3a\Pl + 

H . . d . . . (3) 

- 12d/3+(/?^ _ 3/3!) + 3d-(/?^ + /3!) - |/3+(/3^ - 3/3!)] 

Ri'^R^^ = 3Q01 + /3!)2 + 90a\$l + /?!) + 36d^ + 12a^ 

- 36d(/?+/5+ + /?_/5_) + 601 + i3^_) - 36a{l3l + l3^_) - 3Qaa^ 

Introducing the volume element v = e^", the Einstein field equations yield the equation for 
v: 

?)2 = 3Af 2 + 3Mv + (5) 

Here are A an anisotropy parameter, M the total mass of the dust and A a cosmological 
constant. The equations for P± are 

where a± are constants satisfying + a! = A"^. Thus we can define a new angular variable 
7 by a+ = v4sin(7 — |) and a_ = v4cos(7 — |) that characterizes the different solutions. 

A. The Weyl Curvature Tensor for a classical Bianchi type I model 

In the classical theory we have no factor ordering problems. Inserting the classical 
equations |^ and |^ and expressing them in terms of the volume element v = e^", we get the 
Weyl curvature invariant: 
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3Mv 3Av^ 3Mv 3Av'^ 

2t2\ 1 + — ^ + — ^cos37 + — ^ + 



A^ A^ ' A^ A^ 



(7) 



where the — is for expanding and + is for contracting solutions. Except for degenerate 
cases this curvature invariant diverges as ^ for small values of v. For large values of v 
(if the classical equations allow it) this invariant tends to zero. If we write z = ^ = 



±\/l + ^ + ^ we can write: 



1 fi 

C^^^'C^p,& = ^— (1 - 2ZCOS37 + z") 
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As expected, the Einstein-deSitter case, A = 0, yields a zero Weyl tensor. The Weyl tensor 
is also zero in a few other cases. Let us from now on assume that A ^ 0. Then ifl^ the Weyl 
tensor is zero, then cos 87 = ±1 and z = ±1. IfA = M = these solutions correspond to 
the line element: 

ds"^ = ~dt^ + A^t^dx^ + dy^ + dz^ 
which is easily seen to be locally flat space, if we make the transformation: 

T = t cosh Ax 
X = t sinh Ax 

The last degenerate case where the Weyl tensor is zero, is only a special hypersurface in 
space time. It is the hypersurface where 

M 

(A < 0,M ^ 0). 

All in all we can say that: 

The Weyl curvature invariant, C^^'^^Cafi^s, will for the family of Bianchi Type 
I universes with dust and a cosmological constant, diverge as ^ near the initial 
singularity, except for a set of models of measure zero. 

This seems to come in conflict with the WCC in its strongest form. 

Wainwright and collaborators suggested that it is not the Weyl curvature in- 

variant, C'^^^Co.p^s, that would represent a "Gravitational entropy", but the invariant 

P = — DavD ■ Let us also check out this more restrictive form of the WCC. 
The expression for in a Bianchi type I model is 

2 _ 4 1 + z^ -2z cos 37 
~ 27 M2 + 2AM V + 4A2^;2 

If the universe contains dust, this entity diverges as ^ near the initial singularity. In the 
absence of dust, this entity diverges as the Weyl tensor; the Ricci square is just a constant. 
Even this more restrictive form of the WCC does not seem to agree with the classical Bianchi 
type I universe. 

B. The WCC: Local version 

What conclusions can we draw from these results? Both the entities C°'^'^^Cai3^s and 
diverge near the origin. Thus the local version of the WCC is set in doubt. Based on the 
above investigations we have to conclude: 



^Iff=If and only if 
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If the WCC is ment to be correct in its strongest form, we have to conclude that the 
classical Bianchi type I model, does not give a good description of the initial stages 
of the universe. The models are too simple to take into account the true behaviour of 
the universe. 

The entity suggested by Wainwright and collaborators are more likely to describe a 
"gravitational entropy" . 

Including a positive cosmological constant the universe is unconditionally driven to 
isotropization and — ^ as f — > oo. 



III. THE WEYL CURVATURE CONJECTURE FOR THE LEMAITRE-TOLMAN 

MODELS: LOCAL VERSION 



Choosing universal time gauge we can write the Lemaitre-Tolman (L-T) models as: 

ds"^ = -dt^ + Q^dr^ + R^{de^ + sin^ Odcp'^) (8) 

where Q = Q{r,t) and R = R{r,t). Classically the metric will be a solution of the Einstein 
field equations where the spacetime contains a cosmological constant A and dust. The 
energy-momentum tensor for dust in co-moving coordinates is given by T^,y = diag(p, 0, 0, 0). 
The classical equations will now turn into 

R' = FQ 

where F = F{r) is an arbitrary function, and 



1 



RR^ + -(1-F^)R 



A 
6" 



-R' 



m 



(9) 



The function m = m(r) is given by the integral m(r) = AnpR^R'dr, thus can be inter- 
preted as the total mass of the dust inside the spherical shell of coordinate radius r. The 
square root of the Weyl curvature scalar in the L-T models is: 



2 

7! 



R^^Q^\R~W 



Q'R' 
QR 



R Q QR R^ 
R^ Q~ QR^ R^ 
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Using the results obtained in [|rT| we will investigate the WCC in the L-T models. 



The classical value 



Inserting eq. @ and the classical solution R' = FQ into eq. (|T0D we can write the Weyl 
scalar for all of the L-T models as: 



4^3 



m 
R^ 



m'R 
3R'm 



(11) 
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If we assume that we can write R{r, t) = a(t)r (which is possible in the FRW and the deSitter 
models), the Weyl tensor will be zero for models where 



m oc 



The FRW models do have this behaviour, thus their Weyl tensor vanishes. 

The Schwarzschild black hole case (which is now in Lemaitre coordinates) has m' = 
and the Weyl tensor squared consists solely of the Kretschmann scalar. The square root of 
the Weyl scalar diverges as R~^. 

Let us define a mean dust density function p by the relation: 



4 _ , 

m[r) = -npR (12) 
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Then we can write the Weyl curvature invariant as: 
m{r) is a function of r only, so 

p R 

The physical meaning of the Weyl curvature invariant is now clear; it measures the difference 
between the actual dust density and the mean dust density. \i p = p everywhere, we have a 
FRW universe with a cosmological constant, thus a zero Weyl tensor. Iff the Weyl curvature 
invariant should be ever increasing then the function d = p — p has to obey: dd > 0. 

Let us check out the Weyl tensor more explicitly. Near the initial singularity {R — > 0) 
we can write the solution of the L-T model for all the models approximately as: 

(^rnV {t-to{r))l (13) 
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which leads to 



To avoid intersecting world-lines we have to assume to(r) < [|l^ which makes the difference 
p — p positive. The above entity will therefore diverge as R — > unless = which is the 
FRW case. Thus, the Weyl scalar invariant will diverge when R — > in the generic L-T 
model. 

Let us also check out the scalar suggested by Wainwright and collaborators. The Ricci 
square in the classical L-T model is: 

R'^'R^^ = + 167rpA + 647rV^ 
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We can therefore write 



(15) 



In the small R limit we can use the preceding solution. The A term may be neglected (p 
diverge), so near the initial singularity the entity behaves as 



2mtn 



m'{t - to) 



(16) 



Thus, we have to conclude that (unless tg = 0) 



dt 



< 



(17) 



R — »0 



Bonnor considered the model 1 — F > with a zero cosmological constant upon which 



he studied the behaviour of the same entity. He did however reach the opposite result. But 
as he clearly points out, he inserted the condition Iq = by hand. He assumed an isotropic 
initial singularity, thus demanding that the Weyl tensor is zero initially. As we see, with a 
more general assumption the entity suggested by Wainwright and collaborators will initially 
be a decreasing function of time. 

To see what is actually happening if we push tQ{r) towards zero in the case 1 — > 0, 
we define a new time variable T = t — to{r) . Now a dot means derivative with respect to T. 
Close to T = we can write 

R' = ^mT^g{r,T) 
where g{r,T) ^ 1 — e(r)Ti and e(r) > 0. The entity now becomes 



oc 



{gm') 



2— + mg 



1 2 



(18) 



For to{r) small but non-zero P^ will diverge at the origin. As T increases P^ will rapidly 
decrease, until the g' part becomes dominant in eq. (plsl) . P^ reaches a minimum before 
it will increase again. As tg is driven towards zero, this minimum approaches zero, and at 
tg = 0, the minimum will be at T = where P = 0. We see that the behaviour at the initial 
singularity is drasticly altered near tg = 0. 



Large R limit 

In the large R limit and in the presence of a cosmological constant, we have the deSitter- 
like solutions. As R heads for infinity, we can approximate these solutions as i? oc e^*. The 
ratio ^ is then approximately constant. The Wainwright entity will therefore go as 



P^ oc R-'^ ^ e" 



6Ht 
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It is exponentially decreasing as a function of t, and at the end of the inflationary era, it is 
close to zero. 

In a 1 — = matter dominated universe, we can set A to obtain the solutions of 
eq. (|13D. The behavior of the entity is now exactly the same as for the small R limit, 
but the relation (p!7D will hold for any R. Even in this case the WCC is set in doubt. 



IV. GRAVITATIONAL ENTROPY: REVISITED 



As we have seen the entity diverges badly at the initial singularity, and decreases as 
{t — to)~^ shortly after the initial singularity in the L-T model. Firstly this seems rather 
contra-intuitive. The Weyl tensor squared measures the square of the difference between 
the mean dust density p and the actual dust density p. However, both of these densities 
diverge at the initial singularity. Shortly after the Big Bang both are reduced as the universe 
expands, and p is only reducing as (t — to)~^ compared to p oc (t — to)~^- If the gravitational 
entropy is infinite at the Big Bang singularity, the total entropy also has to be infinite, which 
is contrary to the second law of thermodynamics. 

Let us consider a finite 3-volume V in our space time. According to the first law of 
thermodynamics the matter entropy Sm, the internal energy U and the number of particles 
N in V will evolve as: 

TdSM = dU + pdV - pdN (19) 

If this volume is co-moving then dN = 0. If the matter content is dust then p = 0. There 
are a couple of things to note. Firstly in a dense dust cloud, we expect the internal energy 
of the dust to be high, thus we expect the entropy to be high. Secondly the entropy is an 
increasing function of the volume. Let us therefore consider a co-moving volume V in our 
spacetime. What could the expression for a gravitational entropy be? From the previous 
sections we noticed that an L-T model with p = p is homogeneous. The quantity is an 
inhomogeneity measure in the L-T models. In the absence of a cosmological constant we 
notice that 

The sign of P is here chosen so that the configuration p < p is associated with a P < 



while the more realistic configuration p > p has P > corresponding to [TJ positive. The 
matter entropy is a function of the volume V, thus the gravitational entropy ought to be so 
too. We therefore consider the entity defined by: 

PdV (21) 



V 

Introducing co-moving coordinates we write dV = \/hd^x where is the 3-volume 
element. The integration range is now constant as a function of time and if we integrate 
over a unit coordinate volume which is so small that the integrand is approximately constant, 
we may write 

S = I PdV^PVh 
Jv 
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This volume can in principle be as small as possible and in this sence we can treat the 
entity U = y/hP local entity. This point was emphasized by Hayward Following 



Hayward we can describe entropy by an entropy current vector \1/ with the second law being 

V • ^ > (22) 

where is future-casual. Decomposing ^ into components tangent and orthogonal to the 
material flow vector u {u ■ u = —1): 

^ = su + (23) 

where u ■ = 0, s is the entropy density and is the entropy flux. 

In our case we take the gravitational entropy density proportional to P and assume a 
vanishing entropy flux: 

^ ocPu (24) 
and identify u with the cosmological time vector. Thus 

M-VP + PV-M>0. (25) 
By the identity u ■ Vy/h = \/hV ■ u we get 

M- V(v^P) = M- Vn > (26) 



Thus in this sence we can treat 11 as a local entity. The equation ^ says that the 
directional derivative of 11 along the material flow vector are everywhere positive. In our 
case this simply reduces to 

If we want to relate 11 to the entropy density as above we will expect an increasing 11 as a 
function of time. The difference in "localness" of P and 11 is that P is a scalar or a 0-form, 
while n transforms as the components of a differential form of maximal rank i.e. of a form 
that can be written u = U /\^ dx^ where the index runs over the dimension of the manifold^. 
We will however in this article use the "integrated" version of 11 and therefore continue to 
call the entity 11 non-local even though it the above sense might as well be called local. 

In both articles we considered co-moving dust so we can use these results directly. 

We will therefore investigate the entity 

In the L-T model we can choose the sign ± as mentioned above. We will however consider 
only the case where 11 > 0. 



^Or simply a; = it * P where * is the Hodge star operator. 
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A. Behaviour of 11 in the L-T model 



In the L-T model the entity 11 turns into 

P 



n = ^-R^B! (28) 



In the small R limit we can ignore the A term. If A = 0, 11 is simply 

2i?3 \p\ 



n 



p 



In the limit R — > we get: 



which is finite for t — to = Q and increasing thereafter. This entity has the right behaviour. 

In the deSitter limit with A > the universe is expected to approach homogeneity. In 
this limit the 11 will go as 

m! f p — p 



n = (30) 

Since p will approximately have the same time evolution as p in the deSitter limit, 11 will 
asymptotically evolve towards a constant. The solutions in the deSitter limit can be written 

as i? = exp (ysj^'t + '^f{r)^ where /(r) is some function of r determined by eq. (||). The 
expression for 11 is now 

„ 2 m'^m^ , 

n = fir) 

V3 FA ^ 

The physically realistic cosmologies has f'{r) > thus /(r) is an increasing function of r. 
The relation between /(r), m, F{r) and to(^) is not calculated here but can be found by 
using the result obtained by Zecca |jl5|. The qualitative behaviour is, however, expressed in 
the present relation. 

It appears as if the deSitter state is a stable state. Recall that for a stable thermo- 
dynamical state, the entropy will be constant. Interestingly since 11 7^ in the deSitter 
limit, some of the information of the initial state is present. 

The evolution of the entity 11 for the L-T model can be summarized: 

• Large A and ever expanding: In the initial epoch the dust dominates and 11 is 
increasing linearly in t. The universe is becoming more and more inhomogeneous. 
After the universe has grown considerably, the cosmological constant becomes domi- 
nant, n stops growing and if the cosmological constant is large enough, it decreases 
asymptotically towards a constant value. The universe is smoothened out. 

• Small A and ever expanding: Again the dust dominates initially. The cosmological 
constant is too small to make 11 decreasing. The entity 11 is ever increasing but is 
bounded from above by a relatively large constant value. 
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• Zero A and ever expanding: The 11 will again be ever increasing and will asymp- 
totically move towards a function f{t) = c + bt^ where c and b are constants and p = 3 
iff F2 > 1 and j9 = 1 iff = 1. 

• Recollapsing universe: Due to the dust term, the final singularity will not be 
symmetric to the initial singularity. Hence this entity is asymmetric in time for a 
recollapsing universe. 

In the L-T models we see that 11 behave just like the WCC suggests. 

The Schwarzschild spacetime has m(r) = constant and has a vanishing Ricci tensor. If 
we look at the entity 11 in the region outside the Schwarzschild singularity 11 will diverge. 
This is in some sence the maximal possible value of 11, the Weyl tensor is as large as possible 
and the Ricci tensor is the smallest as possible. Thus at this classical level it seems that 
the Schwarzchild spacetime has the largest possible 11 which is a good thing if one wants to 
connect 11 with the entropy of the gravitiational field. 



B. n in the Bianchi type I model 



We can now use the same entity 11 for the Bianchi type I model which we "derived" for 
the L-T model. The expression for 11 in the Bianchi type I model is 

/- / 1 ^2 _ 2^ cos 37 Y 

~ ~ sT! VM2 + 2AMv + 4A2z;2 ) ^^^^ 

where y/h = v. Near the initial singularity we can make a Taylor expansion of 112 ^^.^^ 
order in v: 

^ 8 , , f fSM 2A\ 



Interpreting m = # as the dust density in coordinate space, we see that as v —>■ 

2A 
3\/3m ' 



n — > 2^ . But n only increases iff 3m2 > 4A. Thus if m 7^ this inequality is satisfied for 



A < 0. Investigating the classical solutions for A > more closely we see that iff 3m2 = 4A, 
the volume element v{t) can be written as a part of the exponential map. In the case where 
3m2 > 4A, v{t) can be written as a part of the hyperbolic cosine and if 3m2 < 4A, v{t) can 
be written as a part of the hyperbolic sine. Thus if the solutions expand too rapidly, the dust 
term does not manage to increase the value of 11 immediately after the initial singularity. 
Also in the Bianchi type I case the deSitter limit will yield a constant value of 11: 

lim n = 



3A5 



In the absence of a cosmological constant the result is simply: 

2A / 9 
n = — ^ (2 + 3mt)(l - cos 37) + -m2t2 
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In this case 11 is monotonically increasing. For large t we have H ^ ^t, which is according 
to the same power-law as in the "flat" and A = case of the L-T model. 

To summarize our investigation of 11 in the Bianchi type I model we can say that the 
entity 11 at the initial singularity is a constant determined by the inverse of the dust density. 
For 3m^ > 4A it will increase immediately after the initial singularity. In most cosmological 
considerations it is assumed that the cosmological constant is small. The exception is in 
the inflationary era in which the vacuum energy dominates over all other matter degrees of 
freedom. The inflationary era will smooth out anisotropics as well as inhomogeneities, and 
the behavior of 11 in this case is therefore expected. In the deSitter limit 11 will asymptotically 
move towards a constant. Large A means small value, while small A corresponds to a large 
value. This is in full agreement with the L-T models. It is also interesting that the entity 
n is very sensitive to different matter conflgurations. This makes it a lot easier to check 
whether the 11 has the right entropic behaviour. 

It also appears to us that the Kasner universe is a sort of "Worst case scenario". It 
has a horrendous initial singularity, a cigar-shaped singularity where the solution space is 
completely disconnected from the isotropic FRW universe. Also the apparently flat case, 
cos 37 = 1, possess a severe and devastating non-smooth singularity at t = 0. In this case 
the singularity is that of a cone and is a true singularity in the smooth structure of the space 
time . The entity 11 is for v = a continuous function in 7 so the cos 87 = 1 case is not 
particularly different from the cases cos 87 7^ 1. It seems that while the Schwarzschild black 
hole case serves as an upper bound in the inhomogeneous case the Kasner universe serves 
as an upper bound for the anisotropic cosmologies. 

In flgure |I| we have plotted 11 as a function of ^ = in the case of a positive cosmological 
constant. By investigating the entity 11 we note that there is one particular case where 11 
is a constant as a function of V. It is the case where 3m? = 4A and cos 87 = | and it 
corresponds to a line element that can be written [0]: 



ds^ = ^ 



2 



8Arj 



2 



+ 



(32) 



In flgure ^ we have plotted 11^ as a function of V in the case of a negative cosmological 
constant. The expanding solutions can be mapped onto the contracting solutions by chang- 
ing the sign of cos 87. Thus if cos 87 = the entity 11 will be symmetric with respect to the 
turning point where the universe turns from an expanding to a contracting phase. We note 
that only if cos 87 > the function 11 will be larger in the contracting phase than in the 
expanding phase. The case where cos 87 = is symmetric with respect to the turning point 
of the evolution of the Bianchi type I universe. Thus one might say that this case is the 
most symmetric of the Bianchi type I universes, while the cases cos 87 = ±1 are the most 
asymmetric cases conserning the time evolution^ 



we look at the spacetimes these solutions correspond to, we note something interesing. The 
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V. QUANTUM COSMOLOGY AND THE WEYL CURVATURE CONJECTURE 



In the previous section we investigated the entity 11 in the classical cosmologies of the 
L-T model and in the Bianchi type I model. We saw that is had a very promising behaviour. 
What can the theory of Quantum Cosmology say about its behaviour? Is a small initial 
value of n more likely than a large value? Again we will use the results from the two papers 



11] and we will start by investigating the L-T models. 



A. The L-T model 



The full wave function for the universe is a linear combination of particular solutions of 
the Wheeler-DeWitt equation: 



where i runs over some index set. In the paper [|Tl[] we considered semi-classical tunneling 
wave functions. The universe was tunneling from a matter-dominated universe which was 
classically confined to a finite size, into a deSitter like universe. After the tunneling across 
the classically forbidden region the universe described vacuum-dominated models, similar 
to an inflationary model of the universe. Thus in the semi-classical approximation the wave 
function will have a superspace current which will flow asymptotically towards deSitter-like 
solutions. 

The actual expectation values from entities like C^^^^Ca/s^s, and 11 are not explicitly 
obtained for these models, because the actual calculations suffer from "endless" expressions 
and highly time-consuming quantities. We will therefore give a more general description of 
the evolution of the Weyl tensor for the L-T models. 

It would be useful to first recapitulate some of the discussion done in []TT|. First of all we 
discussed tunneling wavefunctions in the WKB approximation. In the WKB approximation 
we assume that the wave function has the form '^wkb = exp(±25'), to the lowest order we 
got the Hamilton- Jacobi equation: 



6r) 



pl2. 



2mR-W(l 







(33) 



If we assume that S = J a{r)dr the resulting equation will be the Hamilton- Jacobi equation 
for a point particle with action cr(r) (r is only a parameter). In the Hamilton- Jacobi equation 



case cos 87 = 1 Kasner universe corresponds to the conar-like universe described in §. The (local) 
symmetry group of the spatial section of all these related spacetimes (with matter content) can 
be expanded to M x Sym{E'^) = M x 150(2). This is also the case for the solutions described by 
cos 87 = —1. These solutions lie just on the opposite side of the Kasner circle from the conar-like 
universes cos 87 = 1. Thus these solutions (the conar-like universe cos 87 = 1, and the "anticonar"- 
like universe cos 87 = — 1) have actually a larger (local) symmetry group than the other solutions 
cos 87 7^ ±1. 
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the functional S turn out to be the action at the classical level. Since the classical action 
can be written as an integral over r the assumption S = J a{r)dr is therefore reasonable at 
the lowest order WKB level. We can interpret the action a as the action of a point particle 
moving in a potential V{R) = ^ [—mR + |(1 — F"^)}^ — ^R^^ with zero energy. The 
WKB wavefunction ipwKB for the point particle can then be written i/jwkb — exp(±icr). 
The two WKB wavefunctions can therefore be related by "ifwKB = cxp(J drlnipw^^^). 
Finding first the wavefunction ip we can then relate its WKB approximation with "^wkb 
through "^wKB = exp(/ dr Inipw k b) ■ 

Let us now ask the question: Is it more likely for a universe with small Weyl tensor to 
tunnel through the classical barrier than a universe with a large Weyl tensor? If this was the 
case this would have been in agreement with the WCC. The question is difficult to answer 
in general but we shall make some simple considerations in order to shed some light upon 
it. 

We assume that the dust density near the origin of the coordinates is larger than further 
out. The homogeneous mass function goes as mh{r) — ^irphr^ where ph is a constant. If 
the dust density is larger near the origin of the r-coordinate than for larger values of r then 
m(r) > mh{r) if we demand equality only at r = rmaa;(we only look at closed universes). 
This will not in general change the size of the universe so we can look at the effects from 
m{r) alone. Since m(r) is greater in general for an inhomogeneous universe than for a 
homogeneous universe, we see that the potential barrier will be smaller for an inhomogeneous 
universe than for a homogeneous universe. Thus apparently an inhomogeneous universe will 
tunnel more easily through the classical barrier than the homogeneous universe. Since an 
inhomogeneous universe will have a larger Weyl tensor than an almost homogeneous one, we 
see that universes with large Weyl tensor tunnel more easily than those with a small Weyl 
tensor. 

If we look at the tunnehng amplitude concerning effects from the A term, it is evident 
that larger A will yield a larger tunneling probability. In the initial era inhomogeneities will 
increase the value of 11. We saw that an inhomogeneous state will tunnel more easily through 
the potential barrier than a homogeneous state. The largest tunneling probability amplitude 
thus occurs for universes with a large cosmological constant and large local inhomogeneities. 
Prom a classical point of view the value of 11 initially was large (but increasing thereafter) , 
but as the universe entered the inflationary era the cosmological constant was relatively big 
so the value of 11 at the end of the inflationary era is relatively small. 

In the initial epoch the universe is not believed to be dust dominated. The dust does 
not exert any pressure and dust particles do therefore not interact with each other. A more 
probable matter content is matter which has internal pressure. Even though gravitation 
tends to make the space inhomogeneous, internal pressure from the matter will try to ho- 
mogenize the space. Since dust is the only matter source in our model, the model only 
indicates the tendency for gravitation itself to create inhomogeneities. 

Since the universes tunnel into a deSitter-like state, the cosmological constant will rapidly 
dominate the evolution. The larger the cosmological constant the lower will the entity 11 be 
after the inflationary era ends. 
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B. The Bianchi type I model 



The general solution of the WD-equation for the Bianchi type I models is (no scalar 
field): 

^(v,P)^ J d^k[c(k)^|J^(v)p(k)e'^■^j (34) 

where i/jj^iv) — v'~2Wl,h{2Hv) is a particular solution of the WD equation (with dust), C{k) 
is a "normalizing constant", and p{k) is a distribution function in momentum space. This 
distribution function satisfies the equation: 

d^k\p{k)\^ = 1 

There is still a factoring problem in turning the classical entities to operators. Let us first 
investigate the expectation value of C°'^'^^Ca/3js- Working in the small v limit it does not 

seem possible to avoid the horrendous ^ + which will contribute with a term: 

/^'/^^*^(^ + ^)^ _ 1 J d'k\Cik)\^\^|JJ:iv)\''\pik)\^\k\^ 

~ J d^k\C{k)\^\ijj:{v)\^\p{k)\^ 

Unless wc have a delta-function distribution at A; = 0; p{k) = S'^{k), the contribution from 
this term to the Weyl Curvature invariant will diverge as Thus, we have to conclude 
that, in the small v limit the expectation value of C^^'^'^Cap-yS goes as: 

{C^^^'C^p^s) « ^ (35) 

just as in the classical case. 

Investigating the invariant R^'^R^j,^., we notice that things are not so easy. The Ricci 

I 4 

square also has a term which presumably would contribute with a term. However, looking 
at the classical expression we see that the Ricci square is independent of the anisotropy 
parameter A. This indicates that at the classical level all terms involving the anisotropy 
parameters, have to cancel exactly. This is not the case quantum mechanically. In the 
quantum case operators do not necessarily commute. Hence we may have contributions 
from terms which classically would cancel each other. In other words, the fact that the 
classical vacuum has R^u = 0, does not mean that the quantum vacuum has i^^^ = 0. 

We assume that (^j) is a set of factor-ordering parameters which represents the "true" 
quantum mechanical system in such a way that — represents the classical system. With 
this parameterization of the factor ordering we would expect the Ricci square expectation 
value to be: 
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where fi is some function of v which has the property: f ~ 0, fj{v) ~ constant. Thus for 
small V and C,j the expectation value would behave as 



and 



oc -77777^ ■ constant (36) 



The Weyl square divided by the Ricci squared is in general finite as v — > for a quantum 
system. The actual value at f = is however strongly dependent on the factor-ordering. 
As the factor ordering parameters approach zero, the value will diverge. Quantum effects 
in the early epoch are essential for the behaviour of this entity near the initial singularity. 
As V — > we expect the quantum effect to be considerable, thus expecting that the factor- 
ordering parameters to be large. 

As indicated already in section II the quantum mechanical expectation value of P will be 
lower and presumably finite in the initial singularity. Therefore the expectation value of 11 
is also presumed to be considerably lower in the initial stages than its classical counterpart. 

Comparing different tunneling amplitudes in the Bianchi type I model is difficult and 
more speculative because the Bianchi type I universe has no classically forbidden region for 
A > 0. This causes the lowest order WKB approximation to be purely oscillatory. The 
lowest order WKB wavefunction will therefore be approximately constant. In the paper 
we did however construct under some assumptions a wavefunction which clearly peaked at 
small values of the anisotropy parameter. Thus these wavefunctions predicts universes that 
have a relatively low value of 11. 



C. The effect of factor ordering 

As we saw the results obtained for the early part of the evolution of the universe is 
heavily dependent on the factor ordering. Even if we tried to capture some of the essentials 
of the factor ordering in some parameters, we did not actually know how we should order the 
operators. It appears however that the factor ordering we have used makes our quantities 
more well-behaved near the initial singularity. 



VI. CONCLUSION AND DISCUSSION 

When Penrose first suggested the Weyl curvature conjecture, the conjecture seemed 
reasonable because of the tendency for gravitation to clump matter together and form in- 
homogeneities. Several investigations of the WCC showed that things where not so simple 
and some authors set the WCC in doubt. Most of these authors investigated local entities 
like P^. In a more recent paper by Rothman ^ it appears as if Penrose even tried to with- 
draw the statement. But as we have shown in this paper there is no reason to withdraw 



the statement as long as we stick to non-local entities. In a early work by Husain [ITB] it 



16 



also appears that a non-local entity in a quantum inhomogeneous model, the Gowdy cos- 
mology, is investigated. In that work a vacuum spacetime is considered where he looks at 
the behaviour of the Weyl tensor squared and his results are basically in agreement with 
ours. However entities like 11 and P diverge for inhomogeneous vacuum spacetimes, hence 
for these spacetimes a more careful treatment is needed. Also Rothman and Anninos |]^,|^ 
investigated a non-local entity^ and their conclusion is similar to our^. Even if we cannot 
say that the Weyl tensor is directly linked to the gravitational entropy, we have shown that a 
certain non-local entity which is constructed from it has an entropic behaviour, and reflects 
the tendency of the gravitational field to produce inhomogeneities. 

In the inflationary era the non-local entity 11 asymptotically evolved towards a constant. 
This is believed to be the correct behaviour in an inflationary era. The massive amount of 
vacuum energy tends to smooth out any inhomogeneities or anisotropics which is present 
from the pre-inflationary era. That this entity did not tend to zero (even though the co- 
moving volume expanded exponentially) could be interpreted to tell that 11 contained some of 
the information from the pre-inflationary era. There were some seeds of the inhomogeneities 
left which could be the seeds needed to form galaxies as we see them today. 
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^They used the volume of the phase space, similar to what we do for classical thermodynamical 
systems. 

^When this paper was almost finished a paper by Pelvas and Lake |17] appeared, where they also 
investigated local entities like P^. They also concluded that these entities cannot be a measure of 
gravitational entropy. 
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FIGURES 




V 

(c) COS 37 = 1 

FIG. 1. Plots of n(y) for the Bianchi type I model with different values of cos 87. The cosmo- 
logical constant is set to A = | and in each plot graphs with m = ^ , 1 and 2 are drawn. We have 
also chosen A = |\/3 so that for all of the above cases 11 — > 1 in the limit V — > 00. 
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(a) cos 37 = — 1 



(b) cos 37 = 1 




(c) COS 37 = — I (d) cos 37 = I 

FIG. 2. Plots of n(F)^ for the Bianchi type I model with different values of cos 87. The 
cosmological constant is set to A = — | and in each plot graphs with m = ^ , 1 and 2 are drawn. 
The expanding phase can be mapped onto the contracting phase by changing the sign of cos 87. A 
is chosen to be equal to |\/3. 
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